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Abstract 

Invariant mass spectra of di-electrons stemming from bremsstrahlung processes are calculated 
in a covariant diagrammatical approach for the exclusive reaction Dp — > p sp npe + e~ with de- 
tection of a forward spectator proton, p sp . We employ an effective nucleon-meson theory for 
parameterizing the sub-reaction np — ► npe + e~ and, within the Bethe-Salpeter formalism, derive 
a factorization of the cross section in the form da£>p^p spnpe + e - /dM = da np ^ npe + e - /dMx kine- 
matical factor related solely to the deuteron (M is the e + e~ invariant mass). The effective 
nucleon-meson interactions, including the exchange mesons tt, a, u and p as well as excitation 
and radiative decay of A(1232), have been adjusted to the process pp — > ppe + e~ at energies 
below the vector meson production threshold. At higher energies, contributions from u and p 
meson excitations are analyzed in both, NN and Dp collisions. A relation to two-step models 
is discussed. Subthreshold di-electron production in Dp collisions at low spectator momenta is 
investigated as well. Calculations have been performed for kinematical conditions envisaged for 
forthcoming experiments at HADES. 
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I. INTRODUCTION 



Di-electron production in scattering processes of hadrons at low energies can be de- 
scribed essentially as bremsstrahlung from incoming and outgoing charged particles. Sev- 
eral formulae for bremsstrahlung, obtained within different approximations, have been 
proposed (for a survey of theoretical approaches to bremsstrahlung reactions see e.g. [1]). 
With the focus on intermediate energies, a covariant approach based on an effective meson- 
nucleon theory to calculate the bremsstrahlung of di-electrons from nucleon-nucleon scat- 
tering has recently been presented in [2], continuing and extending the series of previous 
investigations [3]. In this model the effective parameters have been adjusted to describe 
elastic nucleon-nucleon (NN) and inelastic NN — >• NNir processes at intermediate ener- 
gies; besides, the role of excitations of intermediate resonances has been studied within 
this approach and it is found that at intermediate energies the main contribution comes 
from A resonances (see also Ref. [4]), whereas excitations of higher mass resonances can 
be neglected. The role of higher mass and spin nucleon resonances at energies near the 
vector meson (p, uj and 0) production thresholds have been investigated in some detail 
for proton-proton collisions in several papers (see, e.g., Refs. [5, 6] and references therein 
quoted) with the conclusion that at threshold-near energies the inclusion of heavier res- 
onances also leads to good description of data. However, as demonstrated in Refs. [5, 7] 
calculations with a reasonable readjustment of the effective parameters can equally well 
describe the data without higher mass and spin resonances. In contrast, for di-electron 
production in photon and pion induced reactions excitations of low-lying as well as heavier 
resonances can play a role [8] . 

In the present paper we extend the covariant model [2, 7], which is based on an ef- 
fective meson nucleon theory with inclusion of A isobar contributions and vector meson 
dominance, to the exclusive process Dp — > p sp np e + e~ with incoming deuteron D where 
the di-electron e + e~ is detected in coincidence with the fast spectator proton p sp . We 
calculate the cross section of di-electrons produced primarily in bremsstrahlung processes 
which, to some extent, can be considered as background contribution to other, more com- 
plicate processes. In order to preserve the covariance of the approach, the deuteron ground 
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state and the corresponding matrix elements are treated within the Bethe-Salpeter (BS) 
formalism by making use of a realistic solution of the BS equation, obtained within the 
same effective meson nucleon theory [9]. 

Corresponding experiments are planned by the HADES collaboration at the heavy 
ion synchrotron SIS/GSI Darmstadt [10]. The outgoing neutron n and proton p may be 
reconstructed by the missing mass technique. The very motivation of such experiments 
is to pin down the bremsstrahlung component for e + e _ production in the tagged sub- 
reaction np — > np e + e~ [11]. Detailed knowledge of this reaction, together with the 
directly measurable reaction pp — > pp e + e~ , is a necessary prerequisite for understanding 
di-electron emission in heavy-ion collisions. In heavy-ion collisions the di-electron rate is 
determined by the retarded photon self-energy in medium, which in turn is related to in- 
medium propagators of various vector mesons. In such a way the in-medium modifications 
of vector mesons become directly accessible. 

A broader scope is to achieve a refined understanding of the nucleon-nucleon force 
at intermediate energies. Similar to meson production or real bremsstrahlung, the vir- 
tual bremsstrahlung processes probe some off-shell part of the amplitude providing more 
profound insights into the electromagnetic structure of hadrons, e.g., the electromagnetic 
form factors in the time-like region, not accessible in on-mass shell reactions. Another im- 
portant issue of di-electron emission in NN collisions is to supply additional information 
on vector meson production, in particular ui and (f> mesons [7, 12], which is interesting in 
respect to the Okubo-Zweig-Iizuka rule [13] and hidden strangeness in the nucleon. 

Our paper is organized as follows. In section II we consider di-electron production in 
elementary NN collisions by parameterizing the amplitudes by corresponding Feynman 
diagrams. Parameters are adjusted to pp collisions. In section III we extend the approach 
to deuteron-proton collisions by employing the Bethe-Salpeter formalism with a realis- 
tic solution obtained with one-boson-exchange kernel. Within the spectator mechanism 
picture we derive a factorization formula relating the reactions Dp — > p sp pn e + e~ and 
np — > npe + e~. The summary and discussion can be found in section IV. Compilations 
of useful formulae are summarized in Appendices A and B (where a link to two step 
models is outlined), while the Appendix C describes some details needed for deriving the 
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factorization formula. 



II. DI-ELECTRONS FROM NN COLLISIONS 



A. Kinematics and Notation 

We consider the exclusive e + e~ production in NN reactions of the type 

iVi(Pi) + N 2 (P 2 ) - N[(P{) + N&Pfi + e + (h) + e-(k 2 ) (2.1) 

(For an extension towards including hadronic inelasticities in semi-inclusive reactions cf. 
[14].) The invariant eight-fold cross section is 

rfV = /w 1 % A S 1 T(P 1' ^ **' ^ ^' (2 ' 2) 

2yA(s, m 2 ,m 2 ) * sp in S ri - 

where the kinematical factor A is \(x 2 , y 2 , z 2 ) = (x 2 — (y+z) 2 )(x 2 — (y—z) 2 ); the factor \/n\ 
accounts for n identical particles in the final state, |T| 2 denotes the invariant amplitude 
squared. The invariant phase space volume drf is defined as 



* = (p, + p 2 - p; *») n n ^ijtf- (^) 

The 4-momenta of initial (Pi, P 2 ) and final (P{, P 2 ) nucleons are P = (Ep, P) with Ep = 
\Jm 2 + P 2 , an analogous notation is used for the lepton momenta A; 12 ; m denotes the 
nucleon mass, while the electron mass can be neglected for the present kinematics. The 
invariant mass of two particles is hereafter denoted as s with s = (P 1 + P 2 ) 2 ; along with 
this notation for the invariant mass of the virtual photon throughout the paper we also 
use the more familiar notation q 2 with q 2 = s 7 . As seen from (2.3), the cross section 
Eq. (2.2) is determined by eight independent kinematical variables, the actual choice of 
which depends upon the specific goals of the considered problem. In the present paper 
we are mainly interested in studying the invariant mass distribution of the produced 
electrons and positrons. For this sake it is convenient to choose the kinematics with two 
invariants, s 7 = (ki + k 2 ) 2 and s\ 2 = (P{ + P 2 ) 2 , and three solid angles, <iH* , dVt\ 2 and dVL* ± . 
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This corresponds to creation of two intermediate particles with invariant masses ^J~s^ and 
y/sw with their subsequent decay into two final nucleons and two leptons, respectively, 
as depicted in Fig. 1. For each pair, the kinematical variables will be defined in the 
corresponding two-particle center-of-mass (CM) system. This can be achieved, e.g., by 
inserting in Eq. (2.3) the identities 

1 = J ds 7 d 4 P 1 5(P 2 - s 7 ) 5 (4) (P 7 -h- k 2 ), (2.4) 
1 = J ds 12 d 4 P 12 5(P 2 2 - s 12 )5^(P 12 - Pi - Pfi (2.5) 

and rearranging terms in Eq. (2.3) to separate the invariant phase space volumes for the 
"decays" with P 7 = k\ + k 2 and P 12 = P[ + P' 2 , see Fig. 1. With these conventions we 
arrive at 

1 11 1 

d 8 a = — - — -- Y \T\ 2 —ds 12 ds^ (2.6) 

x P 2 (P 7 -> fci + k 2 ) R 2 {P 12 ^P[ + P^) R 2 {P 1 + P 2 -> P 7 + P 12 ), 
where the two-body invariant phase space volume R 2 is defined as 

R 2 (a + b^c + d) = d 4 P c d 4 P d 5^(P a + P b -P c - P d ) 5(P 2 - m 2 c ) 5{P 2 d - m\). (2.7) 



B. Leptonic tensor 

In the lowest order of the electromagnetic coupling (one-photon approximation) the 
di-electron production process is considered as decay of a virtual photon produced in 
strong and electromagnetic NN interactions from different elementary reactions, e.g., 
bremsstrahlung, Dalitz decay, vector meson decay etc. [15]. For such a process the general 
expression for the invariant amplitude squared reads 

\T\ 2 = W uu p^, (2.8) 

where the momentum P 7 of the virtual photon is denoted as q = P 7 = (ki + k 2 ); e is the 
elementary charge. The purely electromagnetic decay vertex of the virtual photon is de- 
termined by the leptonic tensor = J2 3^j v with the current j M = u(ki, si) ^v(k 2 , s 2 ), 

spins 
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where u and v are the corresponding Dirac wave functions for the outgoing electron and 
positron. For unpolarized di-electrons the leptonic tensor reads explicitly 

lp,v = 4 (k lfX k 2v + k lu k 2fX - g^(h ■ k 2 )) . (2.9) 

The electromagnetic hadronic current J M and the hadronic tensor W^ u = X J^Jt ■> be- 

spins 

sides the electromagnetic interaction, also involve the strong interaction between the in- 
teracting nucleons and, consequently, are of a more complicate nature than and l^. 
In virtue of gauge invariance the electromagnetic tensors obey q^l^ = q v l^ u = q u W^ u = 
q u W uu = 0, and one can omit in all terms proportional to q u and q v and write 

= -2(Ak^k lu + Sl g uu ). (2.10) 

It is evident that the leptonic tensor depends solely upon the kinematical variables con- 
nected with the virtual photon vertex (see Fig. 1) and is independent of the variables de- 
termining the nucleon-nucleon interaction. This implies that, due to Lorentz invariance of 
both the amplitude Eq. (2.8) and the corresponding phase space volume R 2 (P-y — > ki + k 2 ), 
one can carry out the integration over the leptonic variables in any system of reference. 
The integration is particularly simple in the CM of the leptonic pair, where q = 0, 



R* = y/H s 7 , n1)/d>s 1 dVL± and all the time like components of vanish. One obtains 



/ 



e 



\„„ 16vre 4 J.J+v 



dniw^—i^ = ^ — . (2.ii) 

± q 3 s 7 K ' 



The remaining integrals can be computed by evaluating each differential volume R 2 also 
in the corresponding CM system. All together we obtain 



da a ern 



2 

dsi 2 dVL* dVL* l2 i 



d Sl ~ 12ss 7 (47r) 5 i 12 7 12 \| s 2 2 A(s,m 2 ,m 2 ) 



A(s, s 7 , s 12 )X(s 12 , m 2 , m 2 ) ^ ^ ^ ^ 



spins 



where <if2* and dVL* 12 are defined in the CM of initial and final nucleons, respectively; a e 
stands for the electromagnetic fine structure constant. 



C. Lagrangians and parameters 

The covariant hadronic current J n is evaluated within a meson-nucleon theory based 
on effective interaction Lagrangians which consist on two parts, describing the strong and 
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electromagnetic interaction. In our approach, the strong interaction among nucleons is 
mediated by four exchange mesons: scalar (a), pseudoscalar-isovector (n), and neutral 
vector (uj) and vector-isovector (p) mesons [2, 7, 16-19]. We adopt the nucleon-nucleon- 



meson (NNM) interaction terms 

Cnn„ = g*NN<P a , (2.13) 

Cnn. = -^NisYd^r&^N, (2.14) 

Cnn p = -g NNp (n^tN$/ - ^Na^rNdT*/} , (2.15) 

C NNu = 9nnu> (ni»N$* - ^-Na^NdT^j , (2.16) 



where N and <P M denote the nucleon and meson fields, respectively, and bold face letters 
stand for isovectors. All couplings with off-mass shell particles are dressed by monopole 
form factors F M = (A 2 M — fi 2 M ) / (A M — k 2 M ), where k 2 M is the 4-momentum of a virtual 
particle with mass fiM- The effective parameters are adjusted to experimental data on 
iViV scattering. At low energies (below the pion threshold) these parameters are rather 
well known and can be taken from iterated T matrix fits of experimentally known elastic 
phase shifts [18]. At intermediate energies, say in the interval 1-3 GeV, it turns out that 
the pure tree level description basing on the above interactions is not able to reproduce 
equally well the energy dependence of the data [2, 20]. Therefore, following [2, 20], we 
take into account an energy dependence of the effective couplings 

5mm -> gNNM(s) = g e~ 1 ^. (2-17) 

In what follows we employ the parameters /, go and Am from [2] which assure a good 
tree level description of the elastic NN -> NN, NN -> NA and inelastic NN -> iViW 
reactions at intermediate energies. Note that problems with double counting of the A 
degrees of freedom are avoided in such an prescription. 

D. Nucleon form factors and gauge invar iance 

The form of the cross section Eq. (2.12) bases essentially on the gauge invariance of 
hadronic and leptonic tensors. This implies that in elaborating models for the reaction 
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(2.1) with effective Lagrangians particular attention must be devoted to the gauge in- 
variance of the computed currents. In our approach, i.e., in the one-boson exchange 
approximation (OBE) for the strong NN interaction and one-photon exchange for the 
electromagnetic production of e + e~, the current J M is determined by diagrams of two 
types: (i) the ones which describe the creation of a virtual photon with q 2 > as pure 
nucleon bremsstrahlung as depicted in Fig. 2 and 3b, c, and (ii) in case of exchange of 
charged mesons, emission of a virtual 7* from internal meson lines, see Fig. 3a. For these 
diagrams the gauge invariance is tightly connected with the two-body Ward-Takahashi 
(WT) identity (see [21-24] and further references therein quoted) 

Q^(p',p) = e ^±lA (<TV) - S-\p)) , (2.18) 

where is the electromagnetic vertex and S(p) is the (full) propagator of the particle. 
It is straightforward to show that, if (2.18) is to be fulfilled, then pairwise two diagrams 
with exchange of neutral mesons and pre-emission and post-emission of 7* (cf. Fig. 2b)) 
cancel each other, hence ensuring q^J^ = 0, i.e., current conservation (see [25]). This is 
also true even after dressing the vertices with phenomenological form factors. However, in 
case of charged meson exchange the WT identity is not any more automatically fulfilled. 
This is because the nucleon momenta are interchanged and, consequently, the "right" and 
"left" internal nucleon propagators are defined for different momenta of the exchanged 
meson. For instance, the contribution to q^J^ from the bremsstrahlung diagrams Fig. 2a) 
and 2d) reads 

u(P[)T NNM S(P 1 - q) IS' 1 (P 1 -q)- S- 1 ^)} u^) % 2 M P' 2 )T NNM u{P x ) + 
uiPDTNNMuiPi) 1 2 s u{P' 2 ) IS' 1 (P 2 ') - S-\P^ + q)\ S (P^ + q) T NNM u(P 2 ) 



u(P[)T nnm u(Pi) 



u(P 2 )T NNM u(P 2 ), (2.19) 



.(^1 A*m) (k 2 Hm) _ 

where h = P' 2 - P 2 , k 2 = P[ - P u S' 1 (P x ) u(Px) = and u^S' 1 (P 2 ') = 0. In fact, 
g M J M 7^ follows. In order to restore the gauge invariance on this level one must consider 
additional diagrams with emission of the virtual photon by the charged meson exchange 
as depicted in Fig. 3a. Then it is easy to show that the contribution from this diagram 



exactly compensates the non-zero part (2.19), and thus gauge invariance is restored. This 
holds for bar NNM vertices without cut-off form factors. Inclusion of additional form 
factors again leads to non-conserved currents. There are several prescriptions of how 
to preserve gauge invariance within effective theories with cut-off form factors [2, 24-26]. 
The main idea of these prescriptions is to include the cut-off form factors into WT identity 
explicitly and to consider the new relations as the WT identity for the full propagators. 
For instance, Refs [25, 26] suggest to interpret the cut-off form factors as an effective 
account of the self-energy corrections and to present the full propagators, entering the 
WT identity (2.18), as the bare ones multiplied from the both ends of the propagator 
line by phenomenological cut-off functions. Formally, all the Feynman rules to calculate 
ladder diagrams as exhibited in Fig. 2 remain unchanged, while in calculation of diagram 
types as depicted in Fig. 3a the meson-nucleon vertex must be multiplied by the square of 
the cut-off form factor Fp INN (k); consequently a factor F^} 1 NN (k) must be included into 
the effective electromagnetic MM7 vertex. In the simplest case the bare mesonic vertex 
rff = (&n* + fan) receives an additional factor [27, 28] 



becoming 



TTi I ^1 Al ^2 A*M /r> on\ 

~ 1 ~ A2 _U2 ~ A2 _U2 y 2 - 2U ) 
ly m ^2 IV M M 



It can be seen that the " renormalized" vertex (2.21) obeys the WT identity for the full, 
renormalized mesonic propagators. In a more general case, on can add to the vertex (2.21) 
any divergenceless term, which obviously does not change the WT identity. Often it is 
convenient to display in the mesonic vertices some terms which assure the WT identity 
and the divergenceless part explicitly, in which case the corresponding vertex reads as 
(see also Ref. [24]) 

q ■ (fa + fa) 



n M = q f 2 (^ i (ki)-A-\ki))+B(k 1 ,h) 



(fan + fan) - <V 



q 2 



,(2.22) 



where A(/c 2 ) denotes the scalar propagator and B(fa, fa) is an arbitrary scalar function. 
In accordance with [26] , the meson propagators are to be multiplied by cut-off form factors 



at both ends of their lines in the diagram Fig. 3a, resulting in 

A(fc2) = Ekmgl. (2.23) 



Note that the above prescriptions for restoration of the gauge invariance in pn collisions 
are valid only if the effective meson-nucleon interaction vertices do not depend on the 
momentum k of the exchanged meson. This is the case for pseudo-scalar ttNN coupling. 
Instead, if the pseudo-vector ttNN coupling (2.14) is chosen then the Fourier transformed 
four divergence of the currents corresponding to diagrams Fig. 2 contains an additional k 
dependence from the derivatives in the irNN vertex (cf. Eq. (2.16)) yielding 

q^jv ~ u(p' 1 )-/ 5 k 1 u(p 1 )A(kl)u(p' 2 )-f 5 k 1 u(p 2 ) - u{p'^ h k 2 u{p x ) A{k\)u{p' 2 )^k 2 u{p 2 ) . (2.24) 

The contribution of the diagram Fig. 3a with the mesonic vertex (2.21) or (2.22) to 
q^J^ reads very similar to (2.24) but in each term both momenta, k\ and k 2 enter, and 
consequently, the gauge invariance can not be completely restored. Within an effective 
meson nucleon theory with interaction Lagrangians depending on derivatives, the gauge 
invariant coupling with photons is introduced by replacing partial derivatives, including 
the NNM vertices, by a gauge covariant form (minimal coupling). Such a procedure 
generates another kind of Feynman diagrams with contact terms, i.e., vertices with four 
lines, known also as Kroll-Rudermann [29] or seagull like diagrams, see Figs. 3b and c. 
We include therefore in our calculations these diagrams and the corresponding interaction 
Lagrangian 

C NNin = - e -^N lb YA,{r^)N (2.25) 
m n 

with 4-potential and charge operator e of the pion. Gauge invariance is henceforth 
ensured. The contact term contributions become particularly important near the kine- 
matical limits. 

All electromagnetic NNj vertices corresponding to the interaction Lagrangian 

£ e N Nl = -e (N^N) A" + ckN (^^) N (2.26) 

with the field strength tensor T VIX = d v A^ — d^A v , and k as the anomalous magnetic 
moment of the nucleon (k = 1.793 for protons and k = —1.913 for neutrons), should 
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also be dressed by form factors. This means that at least two form factors are needed to 
describe electron scattering from on-mass shell nucleons. 

In the more general case of off-mass shell nucleons even eight terms, satisfying the 
necessary symmetry and gauge invariance requirements, with eight scalar form factors 
contribute to the NN'j vertex. At q 2 < and low nucleon virtuality (P 2 ~ m 2 ) it is still 
possible to restrict this set to two effective form factors to describe electron scattering 
from off-mass shell nucleons (e.g., in A(e, e'p) reactions) by modifying (kinematically) the 
iViV7 vertex to satisfy gauge invariance (see for details, [30] and further references therein 
quoted). Unfortunately, for q 2 > information about the electromagnetic form factors 
can be obtained directly only at q 2 > 4m 2 (e.g., from proton-antiproton annihilation into 
an electron-positron pair or the inverse reaction), while the region < q 2 < 4m 2 remains 
unaccessible in an on-mass shell process. This is just the region which includes vector 
meson production and thus, could provide some tests of the validity of the vector meson 
dominance (VMD) [31] model for the electromagnetic coupling to off- mass shell nucleons. 

Various models have been elaborated to calculate the form factors in the time like region 
below the NN threshold (see, e.g., [24, 25, 32]) which basically treat the electromagnetic 
vertex within an effective meson-nucleon theory in terms of photon couplings directly to 
bare nucleons superimposed to the coupling to the meson cloud surrounding the nucleon. 
Besides, one can apply an analytical continuation of form factors based on VMD [31] which 
suggests that the photon first converts into a vector meson which then couples to hadrons. 
This model provides a successful description of the on-mass shell pion form factor. For 
nucleons, VMD predicts a strong resonance behavior of the time like form factors in the 
neighborhood of vector meson pole masses. However, the dipole like behavior of the space 
like nucleon form factor persuades us that VMD is too strong a restriction. In principle, 
the original conjecture of VMD could be augmented by introducing heavier vector mesons 
(see, e.g., [33]) into the parametrization of the nucleon form factor. 

The eight independent form factors can be defined in terms of positive and negative 
energy projection operators as [24] 

T,(P',P)=e £ A 

Pl=±, P2=± 



2m 



FJ? 1 



pi 



A p \ (2.27) 
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where the third form factor F^' p2 is not an independent one, but is connected with -F 1 P 2 1 ' P2 
via the WT identity. The positive (p = +) and negative (p = — ) energy projection 
operators are denoted as A ± respectively. For the half off-mass shell nucleons only four 
terms contribute to the electromagnetic vertex. As mentioned above, in the time like 
region below the NN threshold these form factors are to be computed as loop corrections 
to the bare electromagnetic vertex [24] and/or as analytical continuation of the VMD 
prediction. Microscopical calculations [25] show that the electromagnetic form factors 
are rather sensitive to model assumptions in the region of vector meson pole masses. At 
low energies they depend weakly on q 2 and can be effectively absorbed into the effective 
parameters of the NN^y Lagrangian [2, 34]. 

In the present paper we are primarily interested in di-electron production at inter- 
mediate energies with the mass distributions sufficiently far from the vector meson pole 
masses so that, following [2, 34], we merely put F\ = 1 and instead of F 2 we use the 
anomalous magnetic moment of the corresponding nucleon, i.e., we use the Lagrangian 
(2.26). However, for purely methodological sakes, we also present some results on di- 
electron production at higher energies, where the e + e~ invariant mass covers the p and ui 
pole masses to evidence effects of the present VMD implementation. 

E. A isobar 

Intermediate baryon resonances play an important role in di-electron production in 
NN collisions [2, 4-6, 8, 15, 20, 28, 35]. At intermediate energies the main contribution 
to the cross section stems from the A isobar [2]. Since the isospin of the A is 3/2 only 
the isovector mesons ir and p couple to nucleons and A. The form of the effective AN 
interaction was thoroughly investigated in literature in connection with iViV scattering 
[18, 36, 37], pion photo- and electroproduction [38-42]. The effective Lagrangians of the 
NAM interactions read [36-38]) 




T d a #„N + h.c, 



V Aa T {d?# a p - d a & p } ll3l5 N] + h.c. 



(2.28) 



(2.29) 
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with /anw = 2.13 GeV and /an p = 7.14 GeV [28]. The couplings are dressed by cut-off 
form factors 

2 

(2.30) 



?NAM 



A 2 — P 

Ly NAM K 

where = 1.421 GeV and Atva p = 2.273 GeV [28]. The symbol T stands for the 

isospin transition matrix (see Appendix A), \I/a denotes the field describing the A. Usually 
particles with higher spins (s > 1) are treated within the Rarita-Schwinger formalism in 
accordance with which the A field is a rank-1 tensor (obeying the Klein-Gordon equation) 
each component of which is a 4-spinor satisfying the Dirac equation as well. To reduce 
the number of redundant degrees of freedom the Rarita-Schwinger field satisfies also a 
number of additional subsidiary conditions (cf. [43]). Nevertheless, such a field does not 
uniquely determine the properties of spin-3/2 particles. It is known that an arbitrary 
field of rank 1 provides a basis for a reducible representation of the Lorentz group, which 
can be decomposed into two irreducible representations corresponding to spins s = 1 and 
s = of the vector field. Correspondingly, an arbitrary solution of field equations for ^a 
will be related to spins s = 3/2 and s = 1/2. In order to eliminate the part corresponding 
to s = 1/2 one usually considers spin projection operators [43] acting on an arbitrary 
solution of the field equations which ensure the uniqueness of the description of particles 
with high spins via 

tf« = pf #0 (2.31) 

2 

where ^ is a solution of the spin-| field equations and the spin projection operator is 
defined as 

Pf{p) = - \rfrf - -^J9V - - 7V) • (2.32) 

2 3 3m A 3mA v ' 

Then the propagator for the high-spin particles is constructed in a fully analogous way 
with the case s = 1/2. Remaind that formally the propagator of a Dirac particle with 
s = 1/2 can be expressed as a product of a scalar propagator multiplied by the positive 
energy projection operator A + = p + m evaluated at p 2 7^ m 2 . For the Rarita-Schwinger 
propagator, in order to ensure the propagation of degrees of freedom with s = 3/2, one 
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usually includes also the spin projection operator Ps to obtain 

I) II I A 2 

This propagator is discussed in the literature [44-47] with respect to the fact that for free 
particles the positive energy projection operator A + (p) commutes with the spin projection 
operator P% (P), which is not the case for off-mass shell operators. The adopted order 

2 

of their multiplication is rather a convention than a rule. Another source of ambiguity 
is the convention of what to use in (2.33) as the "particle mass", the on-mass shell value 
?fiA or the off-mass shell invariant mass y/p 1 [34, 44] . Note that different prescriptions for 
the propagator differ by corrections of the order p 2 — m\ which, at intermediate energies, 
could be absorbed in slight readjustments of the effective parameters. Indeed, changing 
the order of the operators in (2.33) we find an almost constant modification of the cross 
sections over a wide range of kinematic variables. 

In our calculations we adopted the prescription of [2, 43], i.e., the propagator is taken 
according to Eq. (2.33). In addition, to take into account the finite life time of A, in 
the denominator of the scalar part of the propagator, the mass is modified by adding 
the width, i.e., ira — * jjia — ^1^/2. For the kinematics considered here the mass of the 
intermediate A can be rather far from its pole value, so that the width, as a function of 
p A , is calculated as a sum of partial widths through the one-pion (A — > Nn) and two-pion 
(A — > Np — > N2-k) decay channels [48]. 

The general form of the ANj coupling satisfying the gauge invariance can be written 
as [34, 39, 49, 50] 

Can, = -^^e^(^)7,75T3^-^^e QM (^)75T 3 (^)^ 

-^^Qa^h.T^d^ + h.c, (2.34) 

e a ,(z) = Qafj, + + -(1 + 4z)A] 7a7M , (2.35) 

where A is a constant reflecting the invariance of the free A Lagrangian with respect to 
point transformations [51]. Since observables must not depend on this parameter, A is 
arbitrary. According to common practice one puts A = — 1. The other parameter, z, is 
also connected with point transformations, however it is a characteristic of the off-mass 
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shell A resonance and remains unconstrained. The meaning of this parameter is that 
every coupling to a spin-3/2 field contains also contributions from couplings with spin- 
1/2 components. Some times z is called the off-mass shell parameter. Investigations of 
the role of the off-mass shell quantity, treated as a free parameter, in different observables 
related to the A show [39, 49] that inclusion of z into the calculations requires a slight 
readjustment of the effective parameters g^ which are also free parameters. This means 
that the parameter z and the effective couplings gi must be simultaneously adjusted to 
given observables. A thorough study of the role of couplings to spin-1/2 particles [34, 38] 
has shown that the dependence on z is rather weak, making the off-mass shell parameter 
redundant (see also discussion in [2]). Basing on this observation, we neglect the off-mass 
shell parameter by merely putting Q afl (z) = g afl . The coupling constants g± are taken as 
in [2], i.e., g x = 5.416, g 2 = 6.612 and g 3 = 7.0. 



F. Results for AT" AT — ► NN 



As mentioned above our effective parameters have been fixed in such a way as to 
reproduce reasonably well the results of the study [2] performed to explain the DLS 
data [52] at low energies. The OBE parameters (listed in Table I) and their energy 
dependence have been taken as in Ref. [2]. Figures 4 and 5 show results of our calculations 
of the mass distribution of di-electrons in pp and pn collisions at two values of the kinetic 
energy, 1.04 GeV and 2.09 GeV, corresponding to those considered in [2]. The dotted 
lines depict the contribution of pure bremsstrahlung processes from nucleon lines, i.e., 
di-electrons are produced solely due to nucleon-nucleon interaction via the one-boson- 
exchange potential. In our actual calculations we include four exchange mesons, n, a, p 
and uj mesons supplemented by a "counter" term simulating a heavy axial vector-isovector 
meson, with the goal to cancel singularities of the pion potential at the origin [2]. 

The dashed lines in Figs. 4 and 5 depict the contributions of the A isobar within 
the same OBE potential. The solid lines represent the total cross section including all 
interferences. It can be seen that for pp collisions almost in the whole kinematical range 
the A contribution dominates. Near the kinematical limits the nucleon contribution 
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becomes comparable with A contributions. At smaller invariant masses, the N and A 
contributions interfere constructively, while near the kinematical limit and at beam energy 
2.09 GeV the interference pattern in pp reactions is just so that the total cross section 
resembles the N or A contributions individually. Another observation is the fact that, 
due to isospin factors the cross section in pn reactions is systematically larger than in pp 
reactions. These results are in agreement with [2]. In pn collisions the inclusion of the 
contact terms amplifies the contribution of pure nucleon diagrams. 

After adjusting the model parameters we proceed and present in Figs. 6 and 7 results 
at energies envisaged in the approved HADES proposal [10] for pp and pn processes. It is 
seen that, except for the absolute values, the behavior of the cross section and the relative 
contributions of A isobars and pure nucleon bremsstrahlung basically does not change 
with energy. However, as seen from these figures, the kinematical range of the di-electron 
mass becomes essentially larger covering also the region of vector meson production, i.e., 
p and u, which have been not yet implemented in the calculations. Consequently, at these 
energies the results presented in Figs. 6 and 7 are to be considered as an estimate of a 
smooth bremsstrahlung background. Effects of p and u> excitations will be considered in 
the next subsection. 

In Fig. 8 the isospin effects in pn and pp reactions are quantified. The dot-dashed line 
illustrates the difference between pp and pn processes in pure nucleon bremsstrahlung, the 
dashed line reflects the isospin effects for the A contribution, while the full line is the ratio 
of the total cross sections. Note the nontrivial invariant mass dependence which prohibits 
the use of simple, constant isospin factors to relate pp — > ppe + e~ and pn — > pne + e~ cross 
sections. 

G. VMD effects 

As seen in Figs. 6 and 7, at beam energies Tki n > 2 GeV the kinematical range of the 
invariant mass M covers the p and ui pole masses (corresponding to 0.768 and 0.783 GeV, 
respectively) so that in this region the vector meson nature of the electromagnetic coupling 
of photons with nucleons can show up. In simple terms, the VMD hypothesis [31] implies 
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that photons couple to hadrons (ir mesons, nucleons etc.) solely via intermediate vector 
mesons, in which case the electromagnetic form factor reads 



2 



W« 2 ) = jjjfp. (2.36) 

Such a behavior of form factors can be obtained in a more rigorous way within an effective 
meson-nucleon theory, like the one used in the present paper. For this purpose one 
should consider additionally effective Lagrangians with electromagnetic couplings of the 
vector mesons (only p and uj in the kinematical region we are interested in) with photons. 
This procedure is not unique and one should pay attention to avoid double counting of 
contributions from the Lagrangian with direct NN'-f coupling. Usually [8, 25, 32, 53] the 
electromagnetic 7p and 70; interaction Lagrangians are added to the NN'-f Lagrangian 
(2.26) and are chosen in the form 

£ pH7 = -1^—^Q P{W) ^ (2-37) 

where Q p<yU) \p is the field strength tensor of the p (u) meson. Note that the Lagrangian 
(2.37) should be considered only together with the Lagrangians (2.16) and (2.26), in which 
case the proton electromagnetic vertex reads 



= ( V - ^ ] + ^ + ^(q 2 )^, (2.38) 



where 



F l (q 2 ) = l + 9NNP 



f, 



pi 



2 f 2 

1 , J fry 9NNuj 1 



M 2 - q 2 g N N P Ml - q 2 



*tf> - 1 + t^mF?' (2 ' 39) 

Note that the vertex function (2.38) obeys the WT identity. The introduced coupling 
constants can be estimated from VMD like pole fits [33, 54, 55] and also from the require- 
ment that at g M — > one has F± — > (1 + r 3 )/2 (see also [24, 53]). As a result, one can 
approximately take 

gNNp _ 1 fpy 9NNc 



fpy 2 f^y !]NNp 
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= 1, (2.40) 



which provides the following form of the form factors 




(2.41) 



(2.42) 



Since both p and uj are not stable the corresponding masses in (2.41) and (2.42) receive 
also imaginary parts, i.e., My — > My — iMyIV/2, where Ty is the total decay width of 
the respective vector meson. In our calculations we take advantage of the fact that a free 
p meson decays mainly into two pions, so that its width, as a function of the invariant 
mass q 2 , can be calculated within the same effective meson-nucleon theory with the result 



has been kept constant T w = 0.03 GeV to simulate the finite resolution [10] (HADES 
envisages an invariant mass resolution AM/M ~ 1% [11]). Note that the VMD model 
can be recovered if, as usually, one takes M p ps M w = My and r p ps T^. 

Some comments are in order here. The VMD hypothesis could be implemented not 
only via the effective Lagrangians (2.37), (2.16) and (2.26) but also by considering the 
simplest form for the coupling [24, 56, 57] 



without the direct term (2.26). The Lagrangian (2.44) corresponds better to the original 
VMD conjecture [31] since it assumes that the electromagnetic coupling NNj occurs 
solely via the vector mesons. It is easily seen that if one employs the Lagrangians (2.16) 
and (2.44), e.g., the form factor F 1 with M p ^ M w = My and T p ps r w , the VMD form 
(2.41) or (2.36) could be obtained [24]. An inclusion of direct terms, like Eq. (2.26), will 
lead to double counting in the corresponding amplitude. Hence, the Lagrangian (2.44) 
simultaneously accounts for vector meson production effects near the corresponding pole 
masses (M p and Mj) and for the background (direct) contribution in the whole kinematical 
range. Consequently a separation of these two kinds of effects is hampered with the VMD 




(2.43) 



where T p (q 2 = M p 



) 0.15 GeV. The width of the to meson in the present calculations 




(2.44) 
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Lagrangian taken as in Eq. (2.44). In our calculations we use the Lagrangians (2.26) and 
(2.37) which allow to distinguish the contribution of direct terms from the vector meson 
production, i.e., the first term in Eqs. (2.41) and (2.42) is referred to as the direct or 
background contribution, while the second one defines the p and u meson contribution 
(see also [25]). 

Figures 9 and 10 illustrate the effects of VMD at two kinetic beam energies, Tj, in = 
2.2 GeV and T^ in = 3.5 GeV, for pp and pn collisions, respectively. In the upper panels 
we present VMD effects for the pure nucleon contribution, while in the lower panels the A 
contribution is included as well. The dashed lines represent the background cross section, 
i.e., the one calculated with only the first term in Eqs. (2.41) and (2.42) (cf. Figs. 6 and 7). 
The dotted lines have been obtained when only the second (p) or the third terms (u) have 
been taken into account. It is seen that the resonance structure of the p contribution is 
not so pronounced being rather broad, because of its relatively low threshold and because 
of the mass dependence of its the decay width. Contrarily, the u contribution has a rather 
sharp, resonance like behavior. The relative contribution of uj and p near the pole mass, 
M ~ M w ~ M p ~ 0.78 GeV , is basically governed by the ratio of meson widths squared, 
r 2 (M )/r 2 (M ) = 0.03 2 /0.15 2 (remind that we attribute to the u meson the actual width 
of 30 MeV to simulate finite detector resolution). 

Here it is worth stressing that within the VMD model the vector mesons are designed 
to mediate the electromagnetic coupling of photons with nucleons. Consequently, they 
contribute in the whole kinematical range of the invariant mass and, except the neighbor- 
hood of the pole masses, are essentially virtual. This implies that, apart from the intervals 
near the pole masses, the cross section can not be presented as a two-step process consist- 
ing of: (i) production of a vector meson resonance with experimentally known width and 
with a mass around the p and/or uj mass, (ii) the independent subsequent decay into a 
di-electron channel (cf. discussion in Ref. [58]). The di-electron emission within the VMD 
model is rather a process of production of a virtual vector particle with the quantum 
numbers of the p or/and uj with subsequent decay into a di-electron (via an intermediate 
conversion into a virtual photon) which contributes in the whole kinematical range of the 
di-electron invariant masses. This issue is discussed in some detail in Appendix B. 
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H. Effects of final state interaction 

Previous studies of threshold-near vector meson production in NN reactions have 
shown [5, 7, 45, 59] that the final state interaction (FSI) between nucleons plays an 
important role. It has been also found that considerable corrections from FSI occur at 
low values of the energy excess, As 5 = y/s — 2m — My, where the relative momentum 
(excitation energy) of the nucleon pair is small. With increasing energy excess the relative 
momentum increases too and FSI effects become less important [5, 7]. 

In reactions of di-electron production the kinematical situation is rather different. As 
seen from Eq. (2.12) the invariant mass, s 7 , of the e + e~ pair varies from the photon point 
to a maximum value dictated by kinematics. Similar to the case of on-mass shell vector 
meson production (cf. Eq. (2.4) in Ref. [7]), in Eq. (2.12) an integration over the excitation 
energy su is to be performed. However, in this case the kinematical range of su and, 
consequently, the range of the relative momentum of the nucleon pair is rather different 
and strongly depends on the value of the di-electron invariant mass. With increasing 
di-electron mass FSI effects are expected to increase. For instance, the results of Ref. [7] 
indicate that for pn reactions at di-electron mass near the u pole mass, FSI effects lead to 
an increase of the cross section by a factor ~ 2, while at lower masses the FSI effects are 
expected to diminish. However, with increasing di-electron mass, the kinematical range 
of S12 shrinks, so that FSI effects are expected to increase. 

To take into account FSI we employ here the same model as in a previous study [7] of 
the vector meson production based on the Jost function formalism [60], which provides a 
good description of NN interaction and phase shifts at low relative momenta. Detailed 
results are displayed in Fig. 11, where FSI effects have been calculated at four values of 
kinetic beam energy. The solid lines illustrate the magnitude of FSI corrections, while the 
dashed lines correspond to results obtained without taking into account FSI (cf. Figs. 7 
and reffiglO). As expected, FSI plays a minor role at low values of the di-electron invariant 
mass but increases at the kinematical limit. For completeness, in the lower row we present 
also results with VMD included, where dotted (dot-dashed) lines denote the cross section 
without VMD, including (excluding) FSI. From this figure we conclude, that throughout 
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the kinematical range of interest, effects of FSI are not too large (20 — 50%), but become 
essential and even dominant at the kinematical limits. The cross section exhibits a similar 
behavior in pp reactions and, therefore, are not exhibited here. 

III. DI-ELECTRONS IN THE PROCESS Dp -► p spP ne+e 

Now we are going to implement the parametrization of the amplitude of the process 
np — > npe + e~ in the exclusive reaction Dp — > p sp pne + e~. The latter process will be 
studied by the HADES collaboration [10, 11] with the above mentioned goals. 

A. Formalism 

Let us consider the reaction 

D{P D ) + p(P0 = p( Psp ) + n(Pi)' + p(P^) + e + {k l ) + e (k 2 ) (3.1) 

within the spectator mechanism. The internal neutron of the deuteron interacts with the 
target proton producing a di-electron as a consequence of bremsstrahlung processes, while 
the detected (forward) proton acts as a spectator. In principle, there could be di-electron 
emission from the detected proton, due to final state interaction effects in the three nucleon 
system. The FSI of the spectator with the active pn system can be estimated within a 
generalized eikonal approximation (see, e.g. [61]) by considering the electro-disintegration 
of the 3 He in processes 3 He(e, e'p)pn. A detailed study of such processes [62] shows that 
in parallel kinematics the FSI effects can be safely disregarded and, consequently, the 
di-electron production off the spectator can be neglected indeed. In the present paper 
the parallel kinematics is guaranteed by the choice of the direction of the detected proton 
in the very forward direction, say at 0° — 5°, as envisaged in the experimental proposal 
[10]. The FSI effects in the active pn pair depend on the di-electron invariant mass. Near 
the kinematical limit the relative momentum of the pn pair becomes small, therefore an 
enhancement of the FSI effects is expected in this region. 

As in the previous section, we choose the kinematics with two intermediate invariant 
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masses s 7 and si 2 , as depicted in Fig. 9. The invariant differential cross section reads 

da = ^^^=^^^= - \T\ 2 — — ds~ ds\2 (3.2) 

2(2tt) 11 v /A(s, M 2 , m 2 ) 6 M ^ pms 2E sp 

x r 2 (p 1 + p n -> p 7 + p 12 ) p 2 (p 7 -> h + k 2 ) r 2 (p 12 -> p; + i*) , 

where is the deuteron mass. Integration over lepton variables can be performed as 
above. The same effective meson nucleon theory as in the previous section is employed to 
parameterize |T| 2 . In order to keep the covariance of the formalism and to use directly 
all the previous results we compute the corresponding hadronic electromagnetic current 
within the Bethe-Salpeter (BS) formalism. The current is now 

J„ = «(Pi, s[)u (P' 2 , s 2 ) O a u (Pi, si) K° (Pn,Ps P ) S-\ Psp )v (p sp , s sp ) , (3.3) 

where u (P, s) and v (P, s) are free Dirac spinors. The operator is a short-hand 
notation for the operators of di-electron production in NN interactions within the adopted 
approach. Actually, represents the set of diagrams depicted in Figs. 2 and 3 with all 
nucleon and photon lines truncated; the BS amplitude for the deuteron with total spin 
projection M. D is denoted as ^^ D , and the (modified) inverse propagator of the spectator 
is 5 -1 (p sp ) = (p sp + m). 

Since our numerical solution for the BS equation has been obtained in the deuteron 
center of mass [9], all further calculations will be performed in this system, i.e., in the 
"anti-laboratory" system. In general, the BS amplitude consists of eight partial compo- 
nents. We take into account here the most important ones, namely the S and D partial 
amplitudes. The other six amplitudes may become important only at high transferred mo- 
menta [63, 64], hence for the present process (3.1) with forward detection of the spectator, 
they may be safely disregarded. Observe that the quantity 

V D = (P n ,p sp ) S-\p sp )v ( Psp , s sp ) (3.4) 

being a four dimensional column in the spinor space, acts as a "deuteron spinor" and 
formally replaces in the deuteron current the corresponding neutron spinor. 
The square of is given by 

I E J » J+ " = l E Tr[(P^ + m)^V D VDO M ]. (3.5) 

M.,spins M., spins 
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A direct evaluation (see Appendix C) yields 

\ E V D V D = 2M D (2n) 3 n D (\ Psp \)j(k n + m), (3.6) 

so that 



6 A 1 

M, spins 



l E JaJ + ^ = 2M D (2n) 3 n D (\p sp \)^Tr[(P^ + m)0^(k n + m)0^ 

M, spins 



= 2M D (2nfn D (|p sp |) - E J" + ) NN > (3-7) 

spins 

where the deuteron momentum distribution is 

^ (IPspl) = ^2 (> (IPspl) 2 + (IPspl) 2 ) (3.8) 

with normalization / no (|p sp |) rf 3 p sp ~ 1. (Note that the normalization is not exactly 
unity because, within the BS formalism, besides the main S and D partial waves other 
components, e.g., the negative P waves, enter in the definition of the total momentum 
distribution. Their contribution is however extremely small [64, 65].) 

Equation (3.7) relates the deuteron cross section with the subprocess of di-electron 
production in proton-neutron interactions via 



oz? da on/ \(s NN ,m 2 ,m 2 )) da np 

2Esp ^7T = 2M D ^ x{sq ^ m1) no (| Psp |) (3.9) 

Remarkable is the factorization of the cross section, for which the needed formulae are 
outlined in Appendix C. 



B. Results 

We have calculated the di-electron production in the exclusive process (3.1) at three 
values of the kinetic energy envisaged at HADES [10], T kin = 1.25, 1.90 and 3.5 A GeV . 
As mentioned above, the effects of the final state interaction of the spectator nucleon with 
the "active nucleon" are minimized within the parallel kinematics, where the spectator 
is detected essentially in the same direction as the incident deuteron with approximately 
the same velocity. In our actual calculations we specify, at each considered energy, three 
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angles for the spectator in the forward direction, 6 = 1°, 2° and 5°. The remaining two 
independent kinematical variables, the momentum of the spectator p sp and the di-electron 
invariant mass M, are considered in a large kinematical range. 

In Fig. 13 the dependence of the cross section on the spectator momentum is exhibited 
for two fixed values of the invariant mass and for 9 = 1°. The three curves (solid, 
dashed and dot-dashed) in each panel correspond to three different beam energies. It is 
seen that at each energy the cross section has a maximum at the spectator momentum 
I Psp I — \ |Pd| ( m the anti-laboratory system this corresponds to |p sp | = 0). The widths 
of the distributions increase with increasing energy, which is merely an effect of the larger 
phase space volume. 

In Fig. 14 the mass distribution is depicted. The three columns (left, middle and right) 
correspond to three different angles (9 = 1°, 2° and 5°), while the three rows (upper, 
middle and lower) specify three kinetic beam energies, Tki n = 1.25, 1.9, and 3.5 AGeV. 
For such kinematical conditions the invariant cross section has been calculated at three 
different values of the spectator momentum, |p sp | = 0.45 \P D \ (dot-dashed curves), |p sp | = 
0.5 |Pd| (solid curves), and |p sp | = 0.55 |Pd| (dashed curves), around the maxima of the 
cross section (cf. Fig. 13). The results in Figs. 13 and 14 do not yet include VMD effects 
in the iVTV subprocess, i.e., the results can be considered as estimates of the background 
contribution. 

Fig. 14 illustrates that the shape of the cross section as a function of the invariant mass 
basically reproduces the one in the pn subprocess. It also seen that at fixed energy and 
angle the effect of variation of the spectator momentum (dot-dashed, solid and dashed lines 
in each panel), apart from decreasing the subprocess's phase space volume with increasing 
momentum, reduces to a factor being proportional to the deuteron momentum distribution 
n D (| Ps P |) (here p sp is the spectator momentum in the anti-laboratory system), as seen 
from Eq. (3.9). At moderate values of the invariant mass, e.g., not too close to the 
kinematical limit, the most favorable conditions for di-electron detection are low angles 
9 ~ 1 — 2° and |p sp | = ^ |Pd|- With 9 increasing the cross section at lower values of 

p sp becomes comparable with the cross section at |p sp | = - \Pd\- This can be explained 



24 



that, in spite of |p sp | < - |Pu|, the deuteron momentum distribution no ( | p sp | ) decreases 
(in anti-laboratory system |p sp | 7^ holds), but the phase space increases faster so that 
the total cross section becomes even larger than at |p sp | = - |Pd| (where in the anti- 
laboratory system |p sp | =0), as the right column in Fig. 14 clearly exhibits. 

Finally, in Fig. 15 we present results with VMD effects implemented at the kinetic 
energy T kin = 3.5 A GeV, for which the kinematical range of the invariant mass covers 
the region of p and uj pole masses. As in the pn subprocess the cross section sharply 
increases near the vector meson poles. The contribution to the peak comes mainly from 
the uj excitation. As can be seen in Figs. 7 and 13 at T^ n = 1.9 A GeV the kinematical 
limit of the di-electron invariant mass is located just in the vicinity of the vector meson 
pole masses. In this region the phase space volume for NN reactions shrinks to zero and 
all possible effects of VMD are masked. However, in the deuteron case the phase space 
volume can be enlarged by considering spectator momenta with velocities smaller than 
the initial one. This implies that in the subsystem of the two active nucleons the total 
energy is larger than in the free NN kinematics. Consequently, effects of sub-threshold 
vector meson production can be observed in this region. The lower the spectator mo- 
mentum the larger kinematical range of the allowed di-electron invariant mass can be 
achieved. However, since the deuteron internal momentum distribution sharply decreases 
with increasing spectator momentum, the subthreshold di-electron production at very 
small I p sp I < I |Pd| ( m anti-laboratory the spectator is backward with increasing |p sp |) 
is prohibited. Therefore, it is clear that at threshold-near energies there should be a 
restricted interval for the spectator momentum within which an experimental investi- 
gation of the subthreshold production of vector mesons can be achieved. At too large 
spectator momenta (|p sp | ~ - |Pd|) the suppression originates from the shrinking phase 
space, whereas too small momenta are restricted by the deuteron's internal momentum 
distribution. 

In Fig. 16 we present our results for the cross section (3.9) at the near-threshold energy 
Tkin = 1-9 A GeV at two values of the spectator momentum detected in the very forward 
direction, 6 = 1°. For such a forward kinematics it is quite easy to estimate the effects 
of enlarging the phase space for the elementary subsystem. At T kin = 1.9 A GeV the 
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deuteron momentum is \Pd\ = 5.36 GeV/c so that, e.g., for the spectator momentum 
|p S p| = 0.25 |Pd|j the momentum of the active neutron before interaction is |p n | ~ 
0.75 |Pd| — 4.02 GeV/c. This corresponds to a kinetic energy of T kin ~ 3.2GeV in the np 
subsystem above the vector meson production threshold, hence the cross section at low 
values of the spectator momenta can leak away into the kinematically forbidden region 
for the free np process, as seen in Fig. 16. The solid lines correspond to the spectator 
momenta |p sp | = 0.25 |Pd| (left panel) and |p sp | = 0.35 |Pd| (right panel) respectively. 
One observes that as far as the invariant mass is not too close to the threshold, the 
cross section for the quasi-free kinematics is much larger than at |p sp | ^ — \Pd\, i-e., the 
suppression caused by the deuteron momentum distribution is more important than the 
effect of enlarging the phase space. In the region close to threshold the quasi-free cross 
section falls rapidly to zero and the subthreshold cross section becomes predominant. It 
is also seen that the contribution of p and uj production (dotted lines) is by one order of 
magnitude above the background. In this case the contribution from the background (the 
first terms in Eqs. (2.41) and (2.42)) can be safely neglected, hence a direct investigation 
of the vector meson production becomes feasible. 

We also have investigated effects of FSI in the reaction Dp — > p sp pne + e~ . As we mainly 
consider the parallel kinematics, FSI of the spectator proton with the active np pair can 
be safely neglected and FSI can be important only in the active np pair. In complete 
agreement with the factorization formula and previous estimates [7, 59] the effects of FSI 
are small at small values of the di-electron mass but become visible and important at 
the kinematical limit. This is illustrated in Fig. 17, where calculations are presented for 
T kin = 1-9 A GeV for two values of the spectator momenta. The solid and dashed lines 
depict results for the spectator momentum |p sp | = | \Pd\, while the dot-dashed and 
dotted lines are for |p sp | = 0.55 |Pd|- The former value of |p sp | corresponds to the quasi 
free kinematics, i.e., to the energy of the active np pair near the vector meson production 
threshold. The latter one determines subthreshold energies even for the np subsystem. 
Thus, for both values of the spectator momenta the vector meson production threshold 
is hardly reached and, consequently, FSI effects here are maximized. At lower values of 
the spectator momentum, as well as at lower di-electron invariant mass, the FSI effects 
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are negligibly small. 



IV. SUMMARY 

In summary we have analyzed different aspects of the di-electron production from the 
bremsstrahlung mechanism at energies envisaged at HADES [10] for the exclusive re- 
actions NN — > NN e + e~ and Dp — > p sp pn e + e~ . To calculate the corresponding cross 
sections we employed an effective meson-nucleon theory with parameters adjusted to elas- 
tic NN and inelastic NN — > NNir [2] reaction data with A isobars included and with 
account of vector meson dominance effects. The performed evaluations of bremsstrahlung 
diagrams can be considered as an estimate of the background contribution, a detailed 
knowledge of which is a necessary prerequisite for understanding di-electron production 
in heavy-ion collisions. Our approach is based on covariant evaluations of the correspond- 
ing tree level Feynman diagrams with implementing phenomenological form factors and 
vector meson dominance effects, with particular attention paid on preserving the gauge 
invariance. The covariance of the approach is achieved by direct relativistic calculations 
of Feynman diagrams for NN collisions and by implementation of the Bethe-Salpeter 
formalism for the Dp reaction. The latter case bases on our previously obtained solution 
of the homogenous Bethe-Salpeter equation with realistic interaction [9]. 

In accordance with previous results [2] our calculations demonstrate that in the region 
of invariant masses far from the vector meson production threshold the main contribution 
to the cross section, in both reactions pp and pn, comes from virtual excitations of A 
isobars. Due to isospin effects the cross section for the reaction pn pn e + e~ is larger 
than the cross section for pp — > pp e + e~ by a factor 1.5 — 3. Note that because of (i) 
contributions of the isoscalar a and u exchange mesons, (ii) differences in the electromag- 
netic coupling in ^p and 772 systems, and (iii) interference effects, the isospin enhancement 
is not ~ 9, as one could naively expect from isospin symmetry considerations. In both 
reactions, pn — > pn e + e~ and pp — > pp e + e~ , the bremsstrahlung cross section exhibits a 
smooth behavior as a function of the di-electron mass. Hence, the bremsstrahlung cross 
section can be considered as background contribution. In the kinematical range close 
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to the vector meson pole masses the cross section has sharp maxima, clearly indicating 
that the di-electron production can be also considered as a tool of testing the validity of 
vector meson dominance. The performed investigations of the vector meson dominance 
effects show that they contribute in the whole kinematical range. Near the pole masses 
our diagrammatical approach provides a form for the cross section resembling a two-step 
model formulae, however it preserves the possibility to trace back the essential differences 
between the two models (see also discussion in [58]). 

Having computed the diagrams for the process pn — > pn e + e~ we implemented them 
for the first time in the reaction Dp — > p sp pn e + e~ with detection of a spectator in the 
forward direction. The cross section for this process has been evaluated in a covariant 
approach based on the Bethe-Salpeter formalism. Relativistic effects are negligible here. 
The Bethe-Salpeter formalism has been used rather for the sake of consistency with the 
covariant diagrammatical approach and for convenience in using the results from calcula- 
tions of the subprocess pn — > pn e + e~ . Within the Bethe-Salpeter formalism we derive a 
factorization formula, i.e., the cross section is cast in a form of a product of two factors, the 
one entirely originating from the deuteron structure and kinematics, the other one being 
exactly the cross section of the subprocess pn — > pn e + e~ . In accordance with the factor- 
ization formula the shape of the cross section reflects the one in the elementary reaction, 
except for some corrections from the deuteron wave function. Apart from these common 
features, an essential difference of the reactions Dp — > p sp pn e + e~ and np pn e + e~ 
can appear. Namely due Fermi motion of nucleons in the deuteron it is possible to find 
such kinematical conditions within those the energy balance is shifted in favor of the el- 
ementary subsystem np, so that subthreshold vector meson production becomes feasible. 
However, the performed analysis shows that the subthreshold cross section is quite low 
due to a strong suppression originating from the deuteron wave function. 

Finally, we found that the effects of final state interaction can be neglected at low 
values of the di-electron mass. With increasing di-electron mass the final state interaction 
effects become more important, in particular at the kinematical limit. 

Our results are presented for kinematical conditions accessible in forthcoming experi- 
ments at HADES [10, 11]. Apart from obtaining valuable information for further under- 
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standing of di-electron emission in heavy-ion collisions, this also offers a tool for inves- 
tigation of the interplay of bremsstrahlung process and vector meson dominance effects, 
i.e., the electromagnetic form factors of nucleons in the time like region not accessible 
in experiments with on-mass shell particles and the properties of the half-off mass shell 
nucleon-nucleon amplitude. 

V. ACKNOWLEDGEMENTS 

Fruitful discussions with H.W. Barz, C. Fuchs, R. Holzmann, J. Knoll and A.I. Titov 
are gratefully appreciated. L.P.K. would like to thank for the warm hospitality in the 
Research Center Rossendorf. This work has been supported by BMBF grant 06DR121, 
GSI-FE and the Heisenberg-Landau program. 

APPENDIX A: ISOSPIN IN THE RARITA-SCHWINGER FORMALISM 

In full analogy with the spin-3/2 space, the wave function for an isospin 3/2-particle, 
e.g., the A, can be written as 

Xfr A =T,( 1X l T N\lrA^e* xX i TN , (AI) 

where = coZ (A ++ , A + , A , A~), Xi rjv * s ^ ne wave function for isospin-1/2, and the 
isospin-1 wave function e A (A = ±1, 0) has the same structure as in Eq. (C5) (see Appendix 
C). 

The definition of the transition matrix T in Eqs. (2.28) and (2.29) reads 

(jr A |T a |ir JV ) = (lAir iV |^r A )(e*) a . (A2) 
Using the completeness of Clebsch-Gordan coefficients one finds 

T+T.3 = 5 a p - ^r a T P . (A3) 
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The needed Pauli matrices for isospin 1 are 

1 -y/l 0^ 

-1 

1 
V% 



T = — 

x Vq 



\ 



l 



T = 

y V6 



\ 



T = 

- 1 z 



(A4) 





1 

1 

Finally, define charge rising (lowering) operators T ±1 = (T x ± iT y )/2 with the following 
properties 



1 
1 

V 00 / 



T +1 \p) ~ |A++>, 
T-i \ P ) ~ |A°>, 



T +1 \n) ~ |A+), 
T_i |n) ~ |A") 



(A5) 



to find 



3 T +> 



TJ~T +1 — — |r + , 



(A6) 



Note that (Tr) = 2 (T +1 r_ + T_ir+) + T z r 2 and 

2 



T+T, 



f 
f 



( 

V 



o 1 
o o 





1 



(A7) 



From these equations one finds that the isospin factor for diagonal (with charge conser- 

2 

vation) or non-diagonal (with rising or lowering of the charge) terms is ±-. Usually 

o 

in order to simplify the notation, one changes the normalization (A3) to 



T+T« = 5, 



a/3 



(A8) 



In this case the factor is included into the coupling constants. In our calculations we 
use the normalization (A8). 



APPENDIX B: TWO-STEP MODEL 



In this appendix we establish a correspondence of our diagrammatical approach to 
models based on a two-step mechanism according to which the cross section for di-electron 
production is expressed as a product of two terms: (i) production of an on-mass shell 
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vector meson with the mass around its pole value and known width, and (ii) the subsequent 
decay via the electromagnetic channel conform known branching ratios. To this end we 
recalculate the cross section Eq. (2.12) within the VMD conjecture and cast it in a form 
close to a two-step model. For definiteness, let us calculate the diagram a) in Fig. 2, where 
now the photon couples the nucleon via an isoscalar vector meson, e.g., the u meson. The 
corresponding part of the amplitude T reads 

T(q 2 ) = u(p' 1 , S ' 1 )d mes (k,P2,S2,P2,s' 2 )S(P 1 -q)T^u(p 1 ,s 1 ) 

x g2 _ M | (j-J -f^u^s^ v(k 2 ,s 2 ), (Bl) 

where formally the operator O mes (k, P 2 , s 2 , P 2 , s 2 ) includes all the exchange mesons, meson 

propagators and contributions from the lower vertex Fig. 2a. For further convenience, in 

the vector propagator in Eq. (Bl) we replace the quantity by (since due to gauge 

Mq q 2 

invariance these terms do not contribute to the amplitude (Bl)). Note that Eq. (Bl) is 
valid at any value of q 2 , which generally, could be quite far from the pole mass Mq. Let 
us introduce a hypothetical on-mass shell particle with the invariant mass sy = s 7 = q 2 
and with quantum numbers as those of uj. Then, for such a particle one can use the 
completeness relation for its polarization vectors i\{q 2 ) (A = 0, ±1) to write 

E £aVKa +M V) = + ^ (B2) 

a « 

and 

Aq 2 ) = E[^(Pi. s i)o mes (A;,p 2 , S2 ,p^ S2 )W-9)(r M -eA + )«(pi ) si) 

A 

1 ( M 2 \ e 2 _ 

x ?^li;j? s( ' ;i ' Sl)fo ' 7) " ( * 2 ' S2) 

1 /a /t 2 \ 2 

= J2^NN^NNv(q 2 ^M 2 ,X) (^.)^u(h, Sl )(Cx-l)v(k 2 ,s 2 ), (B3) 

a q ~ M o \ Jui j q 

where Ann-*nnv (q 2 > ^) corresponds to the production amplitude of a vector par- 

ticle with invariant mass q 2 and polarization A. The obtained formula has almost the 
desired factorized form, however, due the dependence upon A of both terms in (B3), the 
factorization is not yet complete. The exact factorization can be accomplished for the 
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squared amplitude after summation over spins and integrating over the leptonic phase 
space, 



^ E / \T(q 2 )\ 2 R2(q^h + k 2 



spins 



E ZXp'u s'l)dme S (k, P 2 , S 2 , P 2 , S? 2 )S(Pi ~ q) • £ X ) «(Pl, «l) 



spins XX' 



x«( Pl , Sl ) (r< w > • &) 5(Px - g)0+ es (A:, P 2 , s 2 , J* S2 )«(p'i, 



1 


2 e 4 / 


((Z 2 - M 2 ) 


g 4 i 



p 4 /■ / ' M 2 \ 

- / P 2 (g + M $ Uh, k 2 ) (B4) 

Carrying out the integration over P 2 (g — > fci + fc 2 ) and observing that it provides a #aa' 
function one gets 



\ E /|r(? 2 )| 2 i22(?-*i + *2) 



spins 



1 



EE- 4 



NN^NNV 



spins X 

where the quantity 



(gVM 2 ,A 



8 7 r 2 v ^r em ( g 2 ) 



(q 2 - M 2 ) 2 + M 2 T 2 tot (q 2 V 



(B5) 



47raL / M 2 x ~ 



r em (g 2 ) = "" u « J 



(B6) 



plays the role of the electromagnetic decay width of a vector particle into a di-electron 
via intermediate creation of a virtual photon. Then the cross section can be written in 
the form 

1 1 



da 



2(27r) 5 v /A(s,m 2 ,m 2 ) 4 



spins 



X 



•Ann->nnv 



- E / ds 12 R 2 (P 1 + P 2 ^P V + Pi 2 )P 2 (Pi2 -> P{ + P 2 ) 
s^T ern (q 2 )/n 



(^_M 2 ) 2 + M 2 r 2 ot (^' 



(B7) 



where the expression in square brackets can be interpreted as cross section for the creation 



of a vector particle with the mass s 7 = q 2 in a iViV process, 



s 7 r em (g 2 )/7r 



(g 2 -M 2 )2 + M 2 rL(g 2 )' 



(B8) 



It is worth emphasizing that the invariant mass q 2 varies in the whole kinematical range, 
< q 2 < q 2 max . Hence, the form of the cross section (B8) is valid at any initial energy, 
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including deep-subthreshold values, q 2 max <C Mq. Therefore, at such invariant masses the 
cross section can be considered as background contribution to di-electron production. In 
the very vicinity of the pole masses, q 2 — > Mq, the cross section (B8) becomes 



— -<r u {M )- -2 — — (B9) 

ds i (v^-M ) +TL/4 

which exactly coincides with results of the two-step model. Remind that our diagram- 
matical approach can be related to the two-step model only if (i) one considers the cross 
section integrated over the di-electron phase space, (ii) the invariant mass s 7 is not too 
far from the pole masses, (iii) interferences between different kinds of vector mesons, u 
and p, are disregarded, and (iv) the contributions of other diagrams (e.g., with A isobars) 
are neglected. 

APPENDIX C: FACTORIZATION 

The BS amplitudes in the deuteron rest system are of the form [64, 65] 

*M + D ( P mPs P ) = Af(fc n + m)^p-i MM (p sp - m)<Ps(Po, |Pn|), (CI) 

x(p sp - m)0 D (p o , |P n |), 

where k n is an on-mass shell four vector related to the off-mass shell neutron vector P n 
as follows (in anti-laboratory system one has P n = — p sp ) 

k n = (E kn ,P n ), E kn = ^P 2 a + m 2 ; p=(p ,P n ), p = l -M D - E sp . (C2) 

(Ps,d(po, |Pn|) are the partial scalar amplitudes related to the corresponding partial vertices 



as 



|TD ; , G S ,d(P0, Pn , f V 

(PsMPoi \ p n\) = — -5 • (C3) 

Md TP \ 2 
~^sp ~P0 



2 sp / 

In Eq. (CI) M D is the deuteron mass, and the normalization factor is Af 
1 ' 1 



8tt 2E sp (E sp + m)' 
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The components of the polarization vector of a vector particle moving with 4- 
momentum p = (E, p) having the polarization projection M. — ±1, and mass M are 



IM 



;Zm + P- 



v M ' x .\ /(/;„ • .1/) J ' 

where £vi is the polarization vector for the particle at rest with 



(C4) 









f 1 ) 


















" "7! 




' ^"71 


—i 













\ J 




I 1 ; 



(C5) 



The Dirac spinors, normalized as u{p)u{p) = 2m and v{p)v{p) = —2m, read 



/ 



u 



s) = ^Jm + E v 



Xs 

crp 



\ 



( 



Xs 



v(p,s) 



m + Er 



\m + E p J 



m + E, 



-Xs 



V 



Xs 



(C6) 



/ 



where Xs = —i^yXs, and Xs denotes the usual two-dimensional Pauli spinor. Note that 
the denominator in (C3) is zero when one particle (the spectator in our case) is on- 
mass shell. This singularity is apparent since it is exactly compensated by the factors 



(Psp -m) S \p sp ) = p 



2 
sp 



m 



from (CI) and (3.3). In terms of the main BS 



components the "deuteron spinor" Vd (3.4) can be written in as 



Vd = M(pL - m 2 ) <p s Vs - -^<PdV s - -^=<p D Vi 



D 



where 



V s = {kn + m 



E sp + m 




-(E sp + m)(<T£) 
-(P n <r)(<r£) 

-(E sp + m) 3/2 1 



E sp + m 



Xs 



{<Tt)Xs 
(Pn^)(^) 

V E sp + m "-"J 



Xs 



(C7) 



(C8) 
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V, 



D 



3(P n £) I E sp + m -(P n a) \ I 1 
(P n <x) -E sp + m J { 



p 2 

n 



X 



( (ggn) ~ \ 
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/ ip . \3/2 3(P„^) 

(Es P + m) p2 



From (C8) and (C9) one finds 



As 



V E sp + m j 



(P n ff) 

V -(p„<r) o ; 



(C9) 



53 V^ = 3(£ sp + m) 2 (fc n + m), 



53 V^i/ D = 9(£ sp + m) 2 (£;„ + m), 

E Wb = E V^s = -3(£ sp + m) 2 (fc n + m). 



(CIO) 
(Cll) 
(C12) 



Now it is straightforward to obtain (3.6) from Eqs. (CI) - (C12). In Eq. (3.6) the deuteron 
wave functions £/s,d ( I Psp I) are related with the half off-mass shell vertices Gs,d (C3) as 

G,s,d(Po = \M D -E SV , |p sp |) 



U s ,d(\v*p\) = 



(C13) 
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TABLE I: OBE parameters used in the calculations (cf. Eq. (2.17) and Ref. [2]). Note that 

JNNtv = 7, — 9NNtt- 




FIG. 1: Choice of kinematical variables for the process iVi(Pi) + N 2 (P 2 ) -> N{(P{) + N^(P^) + 
e + (ki) + e~(k 2 ). 
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e+ e e+ e 




FIG. 2: Bremsstrahlung diagrams for the process N\ + N2 —> N[ + N' 2 + e + e~ in the one-photon 
and one-boson exchange approximation (thick lines denote propagators of either nucleons or 
baryon resonances). 




a) b) c) 

FIG. 3: Contribution of meson exchange currents (a) and seagull terms (b, c) to the process 
N1 + N2 — > N' l +N' 2 + e + e~ , where N\ and N 2 stand for protons and N2 and N[ denote neutrons. 
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M [GeV] M [GeV] 



FIG. 4: Invariant mass distribution of a e + e~ pair in proton-proton collisions. The dashed 
(dotted) curves depict the contribution of diagrams with bremsstrahlung from jAN (jNN) 
vertices. The solid lines are the results of calculations of the total cross section as coherent sums 
of nucleon and A contributions. The calculations have been performed at DLS energies [52] and 
are directly comparably with results of [2]. 
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FIG. 5: The same as in Fig. 4 but for pn reactions. The dotted curves include also contributions 
from jnNN vertices. 
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FIG. 6: Invariant mass distribution of di-electrons in pp reactions at energies envisaged in 
experiments at HADES [10]. Notation as in Fig. 4. 
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FIG. 7: The same as in Fig. 6 but for pn reactions. 
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FIG. 8: The ratio a(pn —* pne + e~) to a(pp — > ppe + e~) as a function of the invariant mass at 
Tjjin = 2.2 GeV. The dot-dashed and dashed curves exhibit results with only nucleon and A 
contributions, respectively. The ratio of the total cross sections is denoted by the solid line. 
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FIG. 9: Illustration of the role of the VMD effects for the invariant mass distribution of di- 
electrons produced in pp reactions. The left and right panels are for initial energies T^ n = 
2.2 GeV and 

Tkin — 3.5 GeV , respectively. In the upper row results are presented for nucleon 
contributions solely, while in the lower one the total cross sections, including A isobars, are 
shown. The dashed lines present the background contribution (cf. Fig. 6), dotted lines exhibit 
the pure VMD effects, i.e., the separate contribution from p and oo mesons. The solid lines 
illustrate the effects of VMD for the total cross section (bottom panels) and for pure nucleon 
contributions (i.e., without A, top panels). 
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FIG. 10: Same as Fig. 9 but for pn reactions. 
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FIG. 11: Illustration of the role of FSI effects in the invariant mass distribution of di- 
electrons produced in pn reactions. The four panels correspond to the initial beam energies 
Tkin = 1-25 GeV , 1.9 GeV , 2.2 GeV and 3.5 GeV, respectively. The dashed lines present 
the contribution without FSI, while the solid lines illustrate the effects of FSI. In the lower row 
effects of VMD are displayed as well: dot-dashed lines are results of background contribution 
without FSI, dotted lines exhibit the background including FSI. 
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FIG. 13: Dependence of the differential cross section 2E SV> - -= for the reaction Dp —* 

dMa 6 p Bp 

p sp np e + e~ on the spectator momentum at two values of the invariant mass of the lepton pair, 
M = 0.2 GeV (left panel) and M = 0.3 GeV (right panel). The solid, dashed and dot-dashed 
lines correspond to beam energies T^ in = 1.25 AGeV, T^ in = 1.90 AGeV and T^ in = 3.50 AGeV 
respectively. The spectator is assumed to be detected at 6 = 1° in the laboratory frame. 
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FIG. 14: Invariant mass distribution 2E, 



da 



'dMcPpsp 

values of the deuteron beam energy, T^ in = 1.25 A GeV (upper panel), T^ in 



for the reaction Dp — > p sp np e + e for three 

1.90 AGeV 



(middle panel) and T^ in = 3.5 A GeV (lower panel) and three values of the spectator angle 
in the laboratory system, 9 = 1° (left column), 9 = 2° (middle column) and 9 = 5° (right 
column). Dot-dashed, solid and dashed curves correspond to values of the spectator momentum 
|Psp| =0.45|Pd|) I Psp I = 0.50|Pd| and |p sp | =0.55|Pd|, respectively. 
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FIG. 15: Invariant mass distribution 2E ST> -5 for the reaction Dp — > p sp np e + e , at 

dMd 6 p sp 

deuteron beam energy T^in = 3.5 AGeV and three values of the spectator angle in the laboratory 
system, 6 = 1°,2° and 5°. Dot-dashed, solid and dashed curves correspond to values of the 
spectator momentum |p sp | = 0.45 |P_d|, |p sp | = 0.50 |Pd| and |p sp | = 0.55|Pz)|, respectively. 
Effects of VMD have been taken into account. 
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FIG. 16: Effects of sub-threshold vector meson production in the reaction Dp — > p sp pn e + e~ at 
low values of the spectator momentum, |p sp | = 0.25 |P_d| (left panel) and |p sp | = 0.35 |Pd| (right 
panel). Dashed lines depict the background contribution, dotted lines are separate contributions 
from p and uj mesons, and solid lines are for the total cross section. 
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FIG. 17: Illustration of FSI effects for the reaction Dp — ► p sp e + e~ at = 1.9 ^4 GeV 
and two values of the spectator momentum. Solid (dashed) lines correspond to results with 
(without) FSI taken into account for the spectator momentum |p sp | = 0.5 |P_d|, while the dot- 
dashed (dotted) lines correspond to |p sp | = 0.55 |P_d|, respectively. Calculations have been 
performed at three different values of the spectator angle, 9 = 1°, 2° and 5°. In the upper row 
solely contribution from background is displayed, while in the lower row VMD effects have been 
included as well. 
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